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ABSTRACT 

In clusters of galaxies, the reaction of the intracluster medium (ICM) to the motion 
of the co-existing galaxies in the cluster triggers the formation of unique features, which 
trace their position and motion. Galactic wakes, for example, are an apparent result 
of the ICM/galaxy interactions, and they constitute an important tool for deciphering 
the motion of the cluster galaxies. 

In this paper we investigate whether Bondi-Hoyle accretion can create galactic 
wakes by focusing the ICM behind moving galaxies. The solution of the equations 
that describe this physical problem provide us with observable quantities along the 
wake at any time of its lifetime. We also investigate which are the best environmental 
conditions for the detectability of such structures in the X-ray images of clusters of 
galaxies. 

We find that significant Bondi-Hoyle wakes can only be formed in low temperature 
clusters, and that they are more pronounced behind slow-moving, relatively massive 
galaxies. The scale length of these elongated structures is not very large: in the most 
favourable conditions a Bondi-Hoyle wake in a cluster at the redshift of z=0.05 is 
12 arcsec long. However, the wake's X-ray emission is noticeably strong: the X-ray 
flux can reach ~30 times the flux of the surrounding medium. Such features will be 
easily detectable in Chandra 's and XMM- Newton's x-ray images of nearby, relatively 
poor clusters of galaxies. 
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1 INTRODUCTION 

In clusters of galaxies, the interactions of the intracluster 
medium (ICM) with a moving cluster galaxy is expected to 
modify both the local properties of the surrounding medium 
and the galaxy itself. One of the manifestations of such inter- 
actions is the Bondi-Hoyle (B-H) accretion (Bondi & Hoyle 
1944). This physical process can be pictured as follows: as 
the galaxy is moving in the cluster, ICM particles are de- 
flected by the galaxy's gravity and concentrate behind it, 
into the galactic wake. 

Intuitively, one might think that B-H accretion creates 
overdense and cool regions of enhanced x-ray emission be- 
hind the galaxies. As a consequence, the hot interstellar me- 
dia (ISM) of these galaxies look as if they have been disfig- 
ured: instead of being azimuthally symmetric they appear 
elongated, or as if they have a 'plume' of x-ray emission 
attached to them. Such elongated features have now been 
identified in the x-ray images of clusters of galaxies (e.g.. 



around NGC 1404 in the Fornax cluster of galaxies; Jones 
et al. 1997). An up-to-date list of wake candidates can be 
found in Stevens, Acreman & Ponman (1999). 

Bondi-Hoyle accretion is not the only manifestation of 
the ISM/ICM interactions which creates elongated features 
behind moving galaxies. Ram pressure stripping can also 
shape the ISM in such a way that it appears elongated. The 
difference in the nature of the two elongated structures is 
apparent: a B-H wake comprises of ICM material, while a 
ram pressure-induced wake contains galactic material (ISM). 

Unfortunately, instrumentation prior to Chandra and 
XMM-Newton has not generally allowed the separation of 
the pure galactic and the wake components, either spatially 
or spectroscopically. Only the elliptical galaxy M86 has of- 
fered us the opportunity to gain more insight into the nature 
of its wake (Rangarajan 1995). The observational fact that 
the metallicity of M86's wake is higher than the metallicity 
in the surrounding medium, has been used to assign it a 
galactic origin. It seems most probable that in the case of 
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M86, stripping of its ISM is currently in action, and that 
its wake consists mostly of galactic material. No other well 
studied example is currently known. The danger arising from 
the inability to separate the wake's emission from the galaxy 
itself, is that the analysis of the galaxy's X-ray data would 
lead to false conclusions for the characteristics of its X-ray 
halo. The dilution of ISM gas by the ICM gas in the B-H 
wake could, for example, lead to confusing results for the 
metal abundances in elliptical galaxies. 

It should be understood that under certain conditions, 
B-H accretion and ram pressure stripping may occur simul- 
taneously, creating wakes which contain a mixture of both 
galactic and intracluster material. Recent numerical simula- 
tions by Stevens et al. (1999) have shown that this picture 
can be indeed correct. However, it is not clear yet which are 
the effects of each separate process, and how the dominance 
of one process over the other depends on the environmen- 
tal parameters. It is expected, though, that dense ICMs and 
high galactic speeds favour ram pressure stripping. However, 
the results of the B-H accretion cannot be foreseen so easily. 
If we want to understand the action of B-H accretion and 
be able to find the environmental dependences, we have to 
study this process in its first principles. 

The aim of this paper is to disentangle the two physi- 
cal processes by studying the action and results of the B-H 
accretion. We address questions such as: under which condi- 
tions B-H accretion occurs in clusters of galaxies; which clus- 
ters are the best candidates for detecting B-H wakes; and 
whether x-ray observations with the Chandra and XMM- 
Newton observatories can reveal B-H wakes. 

The remainder of this paper is organized as follows: in 
§^ we present the methodology we followed to calculate the 
properties of B-H wakes, section §^ discusses the constraints 
and input parameters imposed by the problem itself. The 
results of the simulations are presented in Finally, in 
§P our results are compared to available observations, and 
simulations of similar physical processes. 



2 CREATION AND EVOLUTION OF A B-H 
WAKE 

Consider a test volume behind a moving galaxy. As the 
galaxy travels in the cluster at the speed of «gai, particles of 
the ICM are deflected by the galactic gravitational poten- 
tial, and directed into this volume. The ICM particles that 
are influenced by the galaxy's attraction, and modify their 
direction of motion, are the ones contained within a cylinder 
of radius equal to the accretion radius (-Race) : 



2GAfgal 



(1) 



(Bondi 1952), where Mgai and Cs are the mass of the galaxy 
and the local speed of sound respectively. 

In time dt, the particles that enter the wake, at a posi- 
tion X along the accretion ax;is, are the ICM's particles which 
were initially in a shell of width db and length ds — v^^idt, 
and had impact parameter h (see Fig 0): 

rf^2acc(a;) = 27rfougainicMd& dt, (2) 

where nicM is the number density of the surrounding 
medium. 




Figure 1. Schematic diagram of the configuration used for the 
calculations of the B-H accretion. The galaxy 'G' is moving in a 
cluster of galaxies through an ICM of number density nicM- 



The particles that start with an impact parameter h 
land in the wake along the accretion axis at a distance x 
from the galaxy. Assuming that the ballistic approximation 
is valid (see §3.1.2) we find that the position x is given by: 



-GA/gai + (G^Mg^ai + ^;|alb')'/' 



(3) 



The velocity, iiin(a;), at the position x on the accretion axis 
relates to the initial velocity of the ICM's particles (ugai) by: 

h 



Vin{x) = Ugal- 



(4) 



The effect of these incoming ICM's particles is i) the 
increase of the internal energy [E\-at{x,t)] of the wake, and 
ii) its confinement by the pressure [Pace (a;, t)]: 



Pacc(x,t) dxdt = ^"^P<(^-;) d\^,,{x) 



(5) 



27r_Rw(x,t) 

where d'^z^^dx) is given by eq. ^ and R^(x,t) is the radius 
of the wake at the position x (Fig. |l|). In eq. v[-^(x,t) is 
the velocity of the incoming particles in the wake's frame of 
reference, and it is given by: 

dR^{x, t) 



v'in{x,t) = Vin{x) + 



dt 



(6) 



As the galaxy travels through the ICM, the constant 
replenishment of the wake with particles and energy causes 
a continuous change to its properties. The evolution of the 
wake is governed by two equations: 

• the conservation of energy: 
dEint{x,t) _ 



dt 



T /• +\ ID („ ,\dVv,{x,t) dEincoui{x,t) 

Lbo\{x,t) - P„(x,t) — 1 — (7) 



dt 



dt dt 
and the conservation of the momentum (p): 

= / [P„{x,t) — PicM — P!,cc{x,t)]da 



(8) 



where a is the surface containing the wake. 

In the above two equations, P^{x,t), and Picm are the 
pressures of the wake and the surrounding ICM respectively. 
In eq. ^, the first term [Lhoi{x,t)] represents the energy radi- 
ated away via thermal bremsstrahlung, and it is calculated 
in the appendix The second term is the work done by 
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the wake to the surrounding medium. The third term is the 
energy added to the wake by the accreted particles : 



dEi„ 



(9) 



where it was assumed that the incoming particles have 
enough time to thermalize with the existing particles in the 
wake, and reach a Maxwellian distribution at the tempera- 
ture of kT^(x,t) (see §^). 



3 SIMULATIONS 

Equations ^ and ^ were solved numerically to predict the 
temperature [kT„{x,t)], and the number density [n^(x,t)] 
of the wake at any position x along the accretion axis. The 
time step of the integration process was constant and such as 
to allow the gas in the wake to reach a Maxwellian equilib- 
rium. Knowing the temperature and density at any time t, 
and position x, the luminosity \Li,oi{x,t)] can be found from 
equation The central surface brightness along the ac- 
cretion axis [Eei-E2(0)], and the wake's su rfac e brightness 

and 



distribution are calculated using eq. A 



A4 respectively. 



The initial conditions for the simulations were chosen 
to comply with the conditions associated with this problem 
(Table |l|) and are derived in the next sections. The chosen 
values for the parameters of the ICM and the galaxies' veloc- 
ity represent the conditions found in poor and moderately 
rich clusters of galaxies, which as will become apparent in 
the next sections, are the fertile environments for B-H wake 
creation. 



3.1 Input constraints and initial conditions 

3.1.1 The mass of the galaxy 

The accretion radius (eq. |^), defines the size of the region 
around the cluster galaxy in which the galactic gravitational 
potential is strong enough to change the direction of motion 
of all the ICM particles which happen to be inside that re- 
gion, and direct them into the wake. Particles approaching 
the galaxy with impact parameters less than i?acc are de- 
flected into the wake, while particles with 6 > -Race continue 
their motion passed the galaxy unaffected. Therefore, an 
apparent, and simple condition for a galaxy to have a B-H 
wake, is that its accretion radius should be larger than its 
size. 

As is obvious from eq. |l], the potential of a wake being 
formed behind a cluster galaxy depends on the velocity and 
the mass of the galaxy. If we assume an average velocity for 
the galaxies in a cluster of Ugai = \/3o", where the velocity 
dispersion o is given by the a — T relation in clusters of 
galaxies: 



(10) 



(White et al. 1997; Wu, Fang & Xu 1998), and a sound speed 
c. of: 



we find that the condition -Race > R^ai gives: 



(11) 
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Figure 2. Galaxy mass-radius constraints for B-H accretion in 
different temperature environments. The lines are for -Race = J?gal 
for 1 (solid line), 2 (dashed line), and 5 keV (dotted line) tem- 
perature clusters. Equation ^ is only valid for galaxies which lie 
in the space between the lines and the y-axis. 
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In Fig. g we plot the permitted ranges of galactic masses 
and radii for a range of temperatures of the ICM. Galaxies 
which lie in the range between the plotted lines and the 
j/-axis will produce large-scale B-H wakes. Figure ^ demon- 
strates that in high temperature environments, only very 
compact objects can have galactic wakes produced by the 
B-H accretion. On the other hand, cool environments are 
more likely to host galactic wakes. 



3.1.2 The velocity of the galaxy 

To derive the condition of eq. |l^, an average galaxy velocity 
was assumed, which corresponds to the velocity dispersion 
a. In reality, the galaxy velocities are distributed around 
this value, so that in a cluster there are galaxies which move 
at lower and larger velocities than a. The accretion radius 
of fast moving galaxies becomes very small, and they are 
not expected to show B-H wakes. It is only the slow moving 
galaxies which can have large-scale wakes. Comparing eq. |l^ 
and |l^, it can be concluded that the subsonic velocity regime 
is the most prolific condition for B-H wakes. 

For the calculations of §^ the ballistic approximation 
was used. The question that arises is whether such an ap- 
proximation is valid in the clusters considered. The ballistic 
approximation assumes that there is no interaction between 
the deflected particles as they stream past the galaxy on 
their way to the wake. This assumption is valid only if the 
kinetic energy of the deflected particles at any distance r 
from the galaxy is larger than their thermal energy. This 
condition translates to: 

'^gal+S —Z-C^. (13) 

r 7 

Substituting in eq. ^ the distance r by eq. |l| we find that 
the ballistic approximation is valid only when: 
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Figure 3. Results of the simulation No. 1 for the two extreme impact parameters (6) used in the simulations: 10 kpc (solid line) and 
Race (dotted line), which correspond to the closest and furthest parts of the wake to the galaxy (0.5 and 15 kpc respectively along the 
accretion axis). The evolution of the temperature, number density, and central surface brightness in the energy range (0.1-10) keV are 
shown in (a), (b), and (c) respectively. 



4i ^ = 0.4c? (14) 



For any galaxy velocity that obeys eq. [ij, the ballistic ap- 
proximation is correct. Clearly, the lowest limit that eq. hj 
defines for t;gai is always lower than the average velocity of 
a cluster galaxy (compare equations |l4| and . 

From the above discussion it is clear that the special 
conditions of the problem indicate that only slow moving 
galaxies have wakes with ICM material. We therefore have 
to restrain the simulations to subsonic regime, which means 
that virialized cluster galaxies should not show the leading 
bow-shock, which is a B-H wake's characteristic when the 
galaxy moves supersonically. We decided to use galaxies ve- 
locities equal to the local speed of sound and 0.632cs, as 
the condition of eq. 14 requires. The additional advantage 
of the subsonic velocities is that no shock waves are formed 
as stated above, which would require a different treatment 
than the one presented here. 



4 RESULTS 

In the following sections we report on the results of the sim- 
ulation runs (see Table |l|). Figure ^ presents the evolution of 
the temperature, number density, and central surface bright- 
ness of the wake for the simulation No. 1. The results are 
shown for a distance of 0.5 kpc (solid line) and 15 kpc (dot- 
ted line) from the centre of the galaxy along the accretion 
axis. 

4.1 The temperature of the wake 

As expected, independently of the environment, the trend in 
temperature is a continuous cooling with time. The change 
of the temperature of the wake at a distance of 5 kpc from 
the galaxy in the three different environments studied is 
shown in fig. ^ 

The signature of the environment is in the rate at which 
the wake is cooling. In lower temperature clusters the wake 
cools down more rapidly than in hotter environments: in a 
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Table 1. B-H accretion : simulations 




Figure 4. Temperature evolution at a distance of 5 kpc from the 
galactic centre in different environments: i) kTiQ]^=l keV (dotted 
line), ii) feTicM=2 keV (dashed line), and iii) fcTicM=3 keV (solid 
line). In all three cases the galaxy velocity was equal to the local 
speed of sound. 



Figure 5. Dependence of the temperature of the wake on the 
velocity of the galaxy, in a cluster of fcTicM=l keV, and njcM = 
1 X lO""^ cm~'^. The temperature is measured at a distance of 
0.5 kpc from the galactic centre, along the accretion axis, and at 
a time of 1 X 10^ yr. 



3 keV cluster, fcTw reaches 10 per cent of the temperature 
of the surrounding medium in 3 x lO'^ yr, while a wake in 
a 1 keV cluster needs approximately half the time to reach 
the same level. This difference is understood because in low 
temperature environments, the external pressure is not as 
large as in hotter clusters: the pressure exerted by the ac- 
creted particles (eq. ^) and Picm are lower in cooler clusters. 
As a result the overdense wake expands more rapidly, and 
cools quicker. 

The temperature of the wake was also found to depend 
on the velocity of the galaxy. Figure ^ shows that the lower 
Vgai, the quicker the wake cools. This finding can again be 
understood in terms of the lower external pressure in cool 
clusters; the lower Ugai is the lower Pace is. 

Along the accretion axis we find a decrease in fcTw. At 
any time, and any environment the temperature at the ex- 
tremes of the wake is lower than in regions closer to the 
galaxy (by 10-30 per cent; see fig. |^. Such temperature 
variations might be measurable with the new x-ray satel- 
lites Chandra and XMM-Newton. 



4.2 The density of the wake 

Not surprisingly, the density of the wake increases con- 
stantly: in any environment the increase is approximately 
one order of magnitude in ~ 10^~^ yr. As fig. |6| shows, in 
higher temperature clusters the density of the wake is lower 
than in cooler ones. Although the number of ICM particles 
accreted into the wake per unit time (d^Zacc) is larger in high 
kTicM environments (because the galaxy velocity is larger 
according to eq. ^ , the number density of the wake in poorer 
environment is larger. This unexpected result can be ex- 
plained because it is the higher flux of particles z^cc / da) 
that is responsible for the high number density in lower tem- 
perature clusters. 

Although the effect is not dramatic, in any environment, 
the wake's number density has its higher value close to the 
galaxy, and decreases with distance from the galaxy along 
the accretion axis. This result can be understood in terms 
of a stronger gravitational potential closer to the galactic 
centre. Finally, Figure ^ shows the dependence of the wake's 
number density on the velocity of the galaxy. 
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Figure 6. Variations of the number density at a distance of 
5 kpc from tlie galactic centre along the accretion axis. The 
results for different temperature environments are shown: i) 
kTiQiji=l keV (dotted line), ii) fcTicM=2 keV (dashed line), and 
iii) feTicM=3 keV (solid line). The galaxy velocity was equal to 
the local sound speed in all three simulations. 
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Figure 8. Comparison of the (0.1 - 10) keV central surface bright- 
ness in different environments: i) A:Ticm=1 keV (dotted line), ii) 
kTiciA='i keV (dashed line), and iii) A;Ticm=3 keV (solid line). 
In all cases the galaxy is moving at the local speed of sound (516, 
730, and 894 km respectively), and the profiles presented in 
this plot correspond to a distance of 5 kpc from the core of the 
galaxy. 
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Figure 7. Dependence of the density of the wake on the velocity 
of the galaxy, in a cluster of fcTjcM=l keV, and njcM = 1 x 
10~^ cm~^. The number density is measured at a distance of 
0.5 kpc from the galactic centre, along the accretion axis, and at 
a time of 1 x lO'^ yr. 



4.3 Surface brightness 

Figure ^ compares the evolution of the surface brightness 
S(o.i-io) kcv(O) in a 1 keV, 2 keV, and 3 keV clusters. 
The way that it was calculated is demonstrated in the Ap- 
pendix Rl. The conversion t o an y other energy ranges can be 
performed by 



applying eq. |A3 
I shows the (0.1- 



As fig. shows the (0.1-10) keV central surface bright- 
ness in any environment increases with time until it reaches 
10-30 times the brightness of the surrounding medium. Af- 
terwards, the wake's emission starts declining until it is im- 



mersed in the background, and the wake becomes unde- 
tectable. The time that wakes 'disappear' from the x-ray 
images depends on the richness of the cluster, with higher 
temperature clusters being able to retain the wake signa- 
tures for longer time. Additionally, regions which are fur- 
ther away from the galactic centre fade away quicker than 
regions closer to the galaxy (fig. ^). The consequence is that 
a wake becomes shorter with the course of time. As can be 
inferred from fig. bl wakes in richer clusters live longer. 



4.4 The length of the wake 

The maximum length that a wake can reach is simply de- 
fined by the ballistic theory as the distance from the galaxy 
(along the accretion axis) which corresponds to an impact 
parameter equal to the accretion radius -Race. However, as it 
was shown earlier, the surface brightness of the wake reaches 
the background level at times that depend on the richness of 
the environment. Therefore, it is apparent that the length 
of the wake depends on the time that it is observed, and 
the parameters of the environment which they are in. Fig- 
ure ^ shows for how long the maximum length of a wake 
is measurable, before the wake is immersed into the back- 
ground emission. Although wakes in richer environments are 
shorter, they are detectable for longer periods of time. 



5 DISCUSSION 

5.1 Comparison with observations 

It has been demonstrated that the Bondi-Hoyle accretion 
alone can give rise to asymmetries in the x-ray images of 
normal galaxies in relatively rich clusters of galaxies. How- 
ever, the scale of these asymmetries is small when compared 
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5.2 Comparison with other simulations 
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Figure 9. The age of the wake when the surface brightness of the 
wake's end becomes equal to the surrounding surface brightness 
versus the maximum length in a 1 keV-cluster (diamonds), 2 keV- 
cluster (triangles) , and 3 keV-cluster (squares) . The filled symbols 
correspond to galaxy velocities equal the local speed of sound 
(simulations No 2, 4, and 6 respectively). The open symbols are 
for lower jjg^i (simulations No 1, 3, and 5). 



to the size of the galaxies: in normal conditions the length of 
a B-H wake cannot exceed ^^20 kpc. In a cluster at a rcdshift 
of z=0.05 this length corresponds to ~ 12 arcsec, just above 
the resolution of the ROSAT HRI detector. This fact clearly 
justifies why wakes have been elusive, and why there arc 
only a few examples reported in the literature. Even with 
the new x-ray satellites Chandra and XMM-Newton, wakes 
will be detectable only in nearby low temperature clusters. 

The best known candidate for a B-H wake might be 
NGC 1404 in the Fornax cluster of galaxies. As seen in 
the ROSAT PSPC image (Jones et al. 1997), its wake 
points away from NGC 1399, which is the central galaxy 
in the cluster. In the PSPC data its length appears to be 
~25 kpc. Prom Rangarajan et al. (1995) we find that at 
the distance of NGC 1404 the temperature and number 
density of the Fornax cluster ICM is fcTicM — 1.1 keV 
and nicM — 0.6 x 10~* cm~^. Simulating a galaxy of 
A/gai = 2.5 X 10^^ Mq moving in an environment defined by 
the properties of the ICM around NGC 1404, we find that 
the length of the wake could be ~17 kpc if the galaxy is mov- 
ing on the plane of the sky at a speed of 300 km s~^, and it 
cannot exceed 20 kpc if the galaxy is moving at 500 km s~^. 
A first comparison of the PSPC data with the results of the 
simulation may suggest that B-H accretion cannot create 
a wake as long as the one observed. However, this incon- 
sistency may indicate that the galaxy is moving at larger 
speeds than the ones assumed here, owing to recent infall 
into the cluster. Such speculations cannot be verified, and 
stringent conclusions cannot be reached before higher spatial 
and spectral resolution data of NGC 1404 become available, 
to accurately constrain the parameters of the wake from the 
x-ray observations. 



There have been several studies of the B-H accretion since 
it was first explored by Bondi & Hoyle (1944) in the context 
of stellar accretion. Hunt (1971) performed quantitative cal- 
culations of subsonic and supersonic accretion flows for the 
case of a point source moving in an adiabatic gas. Sophis- 
ticated 3-dimensional simulations have been performed re- 
cently (RufFert 1996 and references there- in). However, most 
of these simulations were not designed to represent the con- 
ditions encountered in clusters of galaxies. They explored 
either temperature and density regimes which were not rep- 
resentative of clusters of galaxies, or the galaxy velocities 
used were too high. Additionally, the results of these studies 
cannot be easily converted to observablcs, making a compar- 
ison with x-ray observations not straightforward. However, 
these simulations provide a first picture of the results pre- 
sented here. The simulations of a galaxy moving subsonically 
presented by Ruffert (1994) show that a downstream over- 
density is created. The off'sct from the centre of the accrotor 
is small compared to the supersonic examples. Additionally, 
as it is clear from the results of Ruflfert (1994), in the sup- 
sonic case there is no bow-shock in front of the galaxy. On 
the other hand the bow-shock is the most prominent struc- 
ture generated by the B-H accretion onto a supersonically 
moving body. As stated by Ruffert (1994), a B-H wake in the 
subconic case "is the result of the superposition of the slow 
radial flow into the accretor and the relative motion between 
the accretor and the medium". 

Recently, more in depth predictions of the ICM/ISM 
interactions in clusters of galaxies have been presented by 
Stevens et al. (1999). Their simulations allow both dynam- 
ical processes to take place: ICM is deflected and concen- 
trated behind the moving galaxy, and ram pressure strips 
the galactic ISM. The result is the creation of downstream 
density enhancements, which consist of both galactic and 
intracluster material. Unfortunately, this study does not as- 
sess the fraction of each component expected in the wake 
for a variety of environmental conditions. As a result, one 
cannot confidently estimate whether the wakes we observe 
consist of galactic or ICM material. 

However, these simulations give the first indication that 
B-H accretion might dominate in low temperature clusters, 
in agreement with our results, and ram pressure creates 
wakes in richer clusters. If one was to compare the results 
of the present study with the simulations of Stevens et al., 
he would have chosen to use their simulation No. lb, for the 
sake of consistency. Indeed, simulation No. lb presents the 
results of a galaxy/ICM interaction in a cluster with tem- 
perature fcricM=l kcV, and density nicM ~ 4 x 10^* cm~"^. 
These parameters of the ICM correspond well to the con- 
ditions studied here. Generally, there is a good agreement 
between both studies: they both predict an increase of ap- 
proximately two orders of magnitude in the n„, with a si- 
multaneous temperature drop by two orders of magnitude at 
5 kpc along the accretion axis and after 5x 10* yr (see fig. 3 in 
Stevens et al.). The agreement between the two studies may 
lead to the false conclusion that both simulations describe 
the same physical process successfully. However, special care 
must be paid, because the galgixy velocity in Stevens et al. 
is supersonic with Wgai = 960 km s~^. As it was proven in 
the previous sections, such velocity regimes cannot produce 



© 0000 RAS, MNRAS 000, 000-000 



8 Irini Sakelliou 



large-scale B-H wakes, because the accretion radius is too 
small. However, the lack of striking differences make us con- 
clude that B-H dominates in No. lb simulation of Stevens 
et al. 
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6 SUMMARY AND CONCLUSIONS 

The motion of a body through a gaseous medium results in 
the creation of an overdense, cool wake: the gravitational 
attraction of the body deflects the surrounding medium's 
particles, which end up being concentrated behind it. This 
process (Bondi-Hoyle accretion) was first explored by Bondi 
& Hoyle (1944) in the context of stellar accretion. The aim 
of this paper was to investigate if this process is at work in 
clusters of galaxies, which clusters are the best hosts of B-H 
wakes, and if such features can be seen in the x-ray data of 
clusters of galaxies. 

The main results of our calculations can be summarized 
as follows: 

• In clusters of galaxies, the ICM can be concentrated be- 
hind the subsonically moving galaxies, although the results 
of B-H accretion are expected to be more dramatic when 
the galaxies move supersonically. 

• Large-scale B-H wakes can be found only in low tem- 
perature clusters, behind slow moving, and massive galaxies. 

• A stable situation is not reached within reasonable 
timescales: the wake is created, its temperature decreases 
and its density increases constantly. 

• The smaller the temperature of the ICM, the longer the 
wake is. However, wakes in richer clusters live longer and 
are brighter. Such features will be easily detectable with 
the Chandra and XMM-Newton satellites in nearby, poor 
clusters of galaxies. 

• The central, x-ray surface brightness of the wake can 
reach 10-30 times the brightness of the surrounding medium. 

• The properties of a B-H wake depend on the velocity of 
the galaxy. Fast moving galaxies, for example, have hotter 
waJies. On the other hand, galaxies with pronounced B-H 
wakes should not have bow shocks, because their motion is 
subsonic. 



APPENDIX A: BOLOMETRIC LUMINOSITY 
AND SURFACE BRIGHTNESS DISTRIBUTION 

The bolometric luminosity Lboi(a;,i) used in eq. |^ to de- 
scribe the energy radiated away, is found by integrating the 
bremsstrahlung emissivity over the volume of the wake (Vw) 
at any time t: 



Lhoi{x, t) 



1.1/2 

1.8 X 10"'" — 5b / kTl^^{x,t)n^{x,tfdV^ 



(Al) 



If we assume that the emitting region is a cylinder, 
the gas in the wake is uniformally distributed, and it 
emits isotropically, the surface brightness [Eei-E2(^)] at 
any distance 6 from the accretion axis is the integrated 
bremsstrahlung emissivity along a column at the projected 
distance 6 from the accretion axis: 



-1/2, 



r=R^(x.t) 



rdr 



(A2) 



where Iei-E2 the integral which defines the energy range in 
which Eei-E2(^') is calculated, and it is given by: 



^Ei-E2 = / exp 



{-—) 



dv 



(A3) 



Integrating eq. A2 we find that the surface brightness 
at any point 6 is: 



Sei-E2(^) — Sei-E2(0) 



1/2 



(A4) 



R^{x,t) ^ 

where the central surface brightness Eei-E2(0) is given by: 



JE1-E2 



(0) = 13.6 X 10"^**n„r^ 



-1/2- 



7b -fEi-Ea 



R„ix,t) (A5) 
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